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to both beams being polarized. We present a brief numerical evaluation of our 
results by studying the i^-factors and the scale dependence of the NLO cross 
section. 
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1 Introduction 

In the last few years measurements of the spin asymmetries A^ [N = p, n, d) in 
longitudinally polarized deep-inelastic scattering (DIS) have provided much new 
information on the spin structure of the nucleon. Theoretical leading order (LO) 
[1-3] and next-to-leading order (NLO) [1-4] analyses of the data sets demonstrate, 
however, that these are not sufficient to accurately extract the spin-dependent 
quark (Ag = q^ — q^) and gluon {Ag = g^ — g^) densities of the nucleon. This 
is true in particular for Ag{x, Q"^) since it contributes to DIS in LO only via the 
Q^-dependence of gi (or Ai) which could not yet be accurately studied experi- 
mentally. As a result of this, it turns out [1-4] that the x-shape of Ag seems to 
be hardly constrained at all by the DIS data, even though a tendency towards 
a fairly large positive total gluon polarization, J^ Ag{x,Q'^ = 4 GeV'^)dx > 1, 
was found [|^, |^, Q. The measurement of Ag thus remains one of the most in- 
teresting challenges for future spin physics experiments. When selecting suitable 
processes for a determination of Ag, it is crucial to pick those that, unlike gi, 
have a gluonic contribution already at the lowest order. Sticking to polarized 
lepton-nucleon interactions, this implies to consider processes less inclusive than 
DIS. Among those is the production of a (charged) hadron with large transverse 
momentum p^- To obtain a large number of such hadrons, it is expedient to go to 
p/iotoproduction, i.e. to the limit when the (circularly polarized) photon which 
is exchanged between the polarized lepton and the nucleon, is almost on-shell. In 
this way one avoids the suppression of the cross section by the photon propagator. 

As was shown recently [^ , a polarized version of the HERA collider with ^/s ~ 
300 GeV would be a very promising and useful facility for studying polarized 
photoproduction reactions. In particular, two of the conceivable processes, single- 
inclusive hadron production and jet production, show strong sensitivity to the 
polarized gluon distribution of the proton and also appear likely to yield statistics 
good enough for a successful measurement |]5|. In the framework of the LO 
calculation performed in |^, the sensitivity of these reactions to Ag is due in 



the first place to the subprocess 7^ — *> qq, where the arrows denote longitudinal 
polarization. As was stressed in p, and as is well-established in the unpolarized 
case, the (quasi-real) photon will not only interact in a direct ( "point-like" ) way, 
but can also be resolved into its hadronic structure. As far as a determination 
of Ag is concerned, such "resolved" contributions (which appear at the same 
order in perturbation theory as the "direct" piece) are to be considered as a 
background. As was shown in [^], the resolved component is subdominant with 
respect to the direct one in certain regions of rapidity and transverse momentum 
of the produced hadron or jet, thus maintaining the clear-cut sensitivity to Ag 
resulting from the direct piece. Focusing on the other hand on the resolved 
component, the study of polarized photoproduction at HERA might even allow 
a measurement of the parton content of polarized photons in the long run |^ - 
a unique task for HERA which makes the polarization upgrade option of HERA 
appear even more fascinating. 

Polarized photoproduction reactions can also be studied in fixed target experi- 
ments with polarized lepton beam and polarized target, like the future COMPASS 
experiment at CERN, or HERMES at DESY. Among other things, one could look 
for charged tracks with large pt also in these experiments, whereas the energies 
would obviously not be large enough for producing decent jets. The resolved 
component at fixed target energies is expected to be generally negligible. 

In order to make reliable quantitative predictions for a high-energy process 
such as polarized inclusive-hadron photoproduction, it is crucial to extend LO 
studies like the one of ^ to NLO by determining the 0{as) QCD corrections. 
The key issue here is to check the perturbative stability of the process considered, 
i.e. to examine to what extent NLO corrections affect the cross sections and spin 
asymmetries relevant for experimental measurements. Only when the corrections 
are reasonably small and under control can a process that shows good sensitivity 
to, say, Ag at the lowest order, be regarded as a genuine probe of the polarized 
gluon distribution and be reliably used to extract it from future data. The first 



basic ingredient for such an extension to NLO has been provided in the past 
two years by the NLO fits to polarized DIS data mentioned above, which yielded 
spin-dependent nucleon parton distributions evolved to NLO accuracy. Focusing 
on the direct part of inclusive-hadron photoproduction, the calculation of the 
polarized cross section to NLO is then completed by using also (unpolarized) 
NLO fragmentation functions for the produced hadron (as provided in 0), and 
by including the 0{as) corrections to the spin-dependent "direct" subprocess 
cross sections for the inclusive production of a certain parton that fragments into 
the hadron. The calculation of the latter is the purpose of this paper. 

An immediate problem arises here, as the direct part on its own is no longer 
a really well-defined quantity beyond the LO. This is due to the fact that be- 
yond LO collinear singularities appear in the calculation of the subprocess cross 
sections for photon-parton scattering which are to be attributed to a collinear 
splitting of the photon into a qq pair and need to be absorbed into the photon 
structure functions. As the latter only appear in the resolved part of the cross 
section, and since factorizing singularities is never a unique procedure, it follows 
that only the sum of the direct and the resolved pieces is independent of the 
factorization scheme chosen and thus is physical. This has been known for a 
long time from the unpolarized case where the corrections to the direct 0, § and 
to the resolved [§] contributions have all been calculated. Nevertheless, we will 
concentrate in this work only on the corrections to the direct part of the polar- 
ized cross section, mainly because this calculation - albeit already being quite 
involved - is much simpler than the one for the resolved piece. Our results will 
therefore only be the first step in a full calculation of NLO effects to polarized 
inclusive-hadron photoproduction. Despite the fact that they are not complete in 
the sense discussed above, we believe our results to be very important, both phe- 
nomenologically and theoretically: As mentioned earlier, the direct component 
dominates at fixed target energies and also still for the HERA collider situation 
in certain regions of phase space. This means that our NLO results should be 
rather close to the true NLO answer in these cases even if the resolved component 



is only taken into account on a LO basis, which in turn imphes that our NLO 
corrections should already be sufficient to shed light on the question of general 
perturbative stability of the process. We also mention in this context that our 
results for the NLO corrections to the direct hard subprocess cross sections will 
help to obtain or to check those for the resolved ones as the abelian ( "QED-like" ) 
parts of the two are the same. 

The paper is organized as follows: in sec. 2 we present the calculation of the 
0{as) corrections to the direct part of polarized inclusive-parton photoproduc- 
tion. Section 3 is devoted to a brief numerical evaluation of our results for HERA 
and fixed target kinematics. Section 4 contains the conclusions. 

2 Calculation of the NLO corrections to the di- 
rect part of polarized inclusive-parton photo- 
production 

2.1 General framework 

The process we want to study is the single-inclusive production of a hadron h in 
photoproduction in collisions of longitudinally polarized electrons (or muons) and 
protons, i.e. e{Pe)p{Pp) -^ h{ph)X. The NLO expression for the corresponding 
spin-dependent cross section is given by 
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the subscripts "++", "H — " in (|1|) denoting the settings of the hehcities of the 
incoming electron and proton. We have introduced the hadronic variables 

T ~U 

V = l + - , W 
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S = {pe+ Ppf , T = {pe- Phf , U = {pp- Phf , (3) 



and the partonic ones 
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Neglecting all masses, one has the relations 
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where x^ (xp) is the fraction of the longitudinal momentum of the electron (pro- 
ton) taken by the photon (by parton i). Similarly, z is the momentum share that 
hadron h inherits from its parent parton j in the fragmentation process. The 
spin-dependent ( "helicity- weighted" ) parton distributions of electrons and pro- 
tons that appear in the expression (§) for the polarized cross section are defined 
as usual by 

Afrix, M') ^ /J^W(^, M') - f:il\-\x, M') , (6) 

where ft'V\ {x, M^) {f^/i) {x, M^)) denotes the probability at scale M of finding 
parton i with positive helicity and momentum fraction x in an electron or proton 
with positive (negative) helicity. As we only deal with the direct case, the only 
parton type occurring for the polarized electron structure functions is the photon, 
and the A/^ coincide with the spin-dependent Weizsacker- Williams spectrumQ 
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^Yor the resolved contribution, one has A/^ -^ A/| in (0), where A/^ is a convolution of 
A/^ in (R) with the polarized photon structure function for parton type k. 



with rrie being the electron (or muon) mass and Qf^^x ^^^ allowed upper limit on 
the radiated photon's virtuality, to be fixed by the experimental conditions. The 
fragmentation function D'j{z,Mp) in (§), describing the fragmentation process 
j -^ h, is of course the usual unpolarized one since we sum over all polarizations 
in the final state. 

Finally, the spin-dependent LO and NLO cross sections for the subprocesses 
72 -^ jX, dAa!:^^'_^j / dv and dAal^-'_^j/dvdw, which have been stripped of trivial 
factors involving the electromagnetic coupling constant aem and the strong one 
as(/i^), are defined in complete analogy with eq. (|I|). Note that, as indicated 
in (^, dAa ^_^J dvdw will explicitly depend on the renormalization scale /x as a 
result of the renormalization procedure for the NLO virtual corrections, and also 
on the scales M, Mp of the parton distributions and fragmentation functions, 
owing to the factorization of initial and final state coUinear singularities. The 
calculation of the dA.d\^^^-/dvdw is the purpose of this paper. 

To conclude this section, let us note that the expression for the wnpolarized 
cross section for single-inclusive hadron photoproduction is similar to the one in 
eqs. (|l]),(||), taking the sum instead of the difference in (|^) and using unpolar- 
ized subprocess cross sections da^^^^/dv, dal^-'_^j/dvdw and parton distributions 
in (0). The latter correspond to taking the sum instead of the difference in (P), 
and for the electron case the unpolarized Weizsacker- Williams spectrum |T^ is 
obtained from (|^) by replacing 1 — (1 — yY — > 1 + (1 — y)"^. When calculat- 
ing the polarized dAa ^_^J dvdw ^ we will at the same time also determine their 
unpolarized counterparts and compare them to existing analytical results in the 
literature ||]. This will serve as a very good check on our calculation. Further- 
more, the unpolarized cross section is needed when one wants to calculate spin 
asymmetries, defined by 

- E.da^/d^pf, ' ^^' 

which are usually the only quantities directly accessible to experiment. 



2.2 LO contributions 

The subprocesses contributing to Ao-^jl^^- are 

7^ ^ ?(g) • (9) 

Here it is understood that the final-state particle in brackets is unobserved and 
integrated over its entire phase space, while the other fragments into the hadron. 
Note that the last process in (P) is symmetric under exchange of q, q. The 
corresponding spin-dependent cross sections read: 
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where Cp = 4/3, Tr = 1/2, and e^ is the fractional charge of the quark. 

2.3 NLO contributions 

Apart from the generic inclusive processes •yq ^ g, 'jq ^ q, and •yg ^ q that are 
already present at the LO level, there are also contributions that can arise only 
beyond the Born approximation. These are 'yg ^ g, 'jq ^ q, and 7g — ^ q', where 
in the latter process q' denotes a quark (or antiquark) of fiavour different from q. 
This means that the following explicit subprocesses have to be evaluated: 

a) the interference between the Born graphs ^q ^ gil): IQ ^ l{g)^ ig ~^ li.Q) 

and the virtual corrections to them, 

b) the real corrections to the Born graphs, 7g — ^ g{(ig), ig ^ li^ig), and 

(gg) 

^q-^q{ {qq) 
{q'q') 



(note that for the latter contribution a finite answer is obtained only if all 
three subprocesses are added), 

c) 19-^9{qq), iq -^ q{qq), iq ^ q'{qq') ■ 

2.4 Regular izat ion of singularities 

It is well known that singularities are encountered when calculating the loop 
diagrams or when performing the phase space integrations for the unobserved 
partons in the 2^3 processes: first of all, the loop-diagrams contain ultraviolet 
divergencies which are removed by renormalization. Adding the renormalized 
loop and the corresponding 2^3 contributions, the infrared singularities which 
are individually present in both parts, also cancel out, and one is left with coUinear 
singularities which are finally removed by the factorization procedure (for the 
contributions from c) only singularities of the latter type occur). Of course, for 
being able to handle the singularities, one has to choose a consistent method 
of regularization. In our calculation we use dimensional regularization for this 
purpose, where (i = 4 — 2e, which is the most convenient and customary choice. 

The calculation of the spin-dependent squared matrix elements requires pro- 
jection onto definite helicity states of the incoming particles (which are taken to 
have momenta pi, P2); which is achieved by using the relations 

m(Pi, hg)u{pi, hq) = - ^1(1 - /ig75) (11) 

for incoming quarks with helicity hg (analogously for antiquarks) and 



e'^ip2,Xg)e*^iP2,X,) ^ 



2(1 -e) 



Pi ■P2 



2pi • P2 



+i^z-^^'%y2Pi 



(12) 
for incoming gluons with helicity A^. The parts independent of hq and A^ con- 
tribute to the unpolarized matrix elements, for which the averaging of gluon spins 
in d dimensions should be performed by dividing by the d — 2 = 2(1 — e) possible 
spin orientations, as has been made explicit in eq. (|T^). 
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As is well known, the use of 75 and the Levi-Civita tensor appearing in 
(|lT]),(|12D is not entirely straightforward in d 7^ 4 dimensions. For our calculations 
we will use the original prescription of 't Hooft and Veltman [|ll|, afterwards sys- 
tematized by Breitenlohner and Maison []T2[ (HVBM scheme), which is usually 
regarded as the most reliable scheme in the sense that its internal algebraic con- 
sistency is well established. In this scheme explicit definitions for 75 and e^^p^ 
are given. In particular, 75 = 'ie^'^''°"7^7i,7p7o-/4!, the e-tensor being regarded as a 
genuinely four-dimensional object with its components vanishing in all unphysical 
dimensions. In this way the (i-dimensional Minkowski space is explicitly divided 
into two subspaces, a four- dimensional one and a (d — 4)-dimensional one, each 
of them equipped with its metric tensor. As a result, apart from (i-dimensional 
scalar products p- q (the usual Mandelstam variables), also their respective "sub- 
space" counterparts can show up in calculations, which renders the calculation 
of traces and phase space integrations somewhat more complicated. Fortunately, 
we can rely in our calculation to a certain extent on known results, as will be 
discussed in the next subsection. 

2.5 Virtual corrections and 2^3 matrix elements 

In |]13[ the NLO corrections to the ( "non- fragmentation" ) part of the hadronic 



single-spin cross section for the production of circularly polarized prompt photons, 
i.e. the QCD corrections for pp — > 7X, were calculated. This calculation involved 
the virtual corrections to the Born graphs qg -^ 7g, gq -^ 7g, qq -^ 7(7, as well 



as the 2 — > 3 matrix elements ab —* jcd. These ingredients were obtained in |]T3 
in the HVBM scheme. We therefore can get the virtual corrections for ^q —>■ gq, 
7g -^ qg, 7^ — > qq and the 2^3 cross sections 7a -^ bed by appropriately 
crossing the polarized photon with the unpolarized incoming parton in the results 



of 13|, which greatly facilitates the calculation. 



The virtual corrections obtained in this way read in the MS scheme: 
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where iVc = 3, 6o = llA^c/6 — nf/3 {uf being the number of active flavours) 
and n is the renormalization scale. Furthermore, 



5T 
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il-v^)/v 



STg^ = —vi/v — v/vi 



(16) 
(17) 



with f 1 = 1 — f . Note that the result for q'-f -^ qg can also be obtained from the 



one of |T3] for qg -^ ^q after crossing and correct adjustment of colour. 



The integration of the real 2^3 matrix elements over the phase space of 
the unobserved particles has been discussed in detail in |jT5|, |14| and needs not 
be recalled here. The technical complications related to the use of the HVBM 
scheme discussed above have been solved in HT^ 
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Adding the renorinalized virtual and the real contributions, all infrared singu- 
larities cancel out. In the next section we briefly recall the factorization procedure 
which removes the remaining coUinear singularities. 

2.6 Factorization 



The factorization procedure based on the factorization theorem [T^ has been 



outlined in, for instance, refs. |T3|, [I^. The mass singularities arise when either 
an incoming particle coUinearly emits another particle (or splits into a pair of 
coUinear particles), or when the "observed" final state particle is collinear to an 
unobserved one. The singular terms attached to the initial legs are separated off at 
the factorization scale M and absorbed into the initial-state parton distributions 
which then obey NLO QCD evolution equations. In particular, if the singularity 
results from a collinear splitting 7 — > gg, it is absorbed into the "pointlike" 
part of the photon structure function. Of course there is freedom in choosing 
the factorization prescription, i.e. in subtracting finite pieces along with the 
pole terms. As already pointed out in the introduction, this is the reason why a 
separation of direct and resolved contributions to a photoproduction cross section 
becomes, strictly speaking, meaningless beyond LO. Final state singularities are 
factorized at the scale Mp into the (NLO) unpolarized fragmentation functions 

T-)h 

As an example, let us briefly discuss the factorization of the polarized q^ -^ 
g{qg) subprocess. This is performed in the easiest way by adding a "counter cross 
section" []T3| which, taking into consideration all possible collinear configurations. 



has the form 

rfAo-AAo a. 
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Uo 
where the dAa'^^~*'^'^{s, v)/dv are the polarized d- dimensional 2 -^ 2 cross sections 



for the processes ah -^ cd, to be found for the HVBM scheme in [0]. Furthermore, 

(A)i7,,(z, M') ^ -i(A)P,,(^) (^)' + (A)/., W , (19) 



where 1/e = 1/e — 7^; + ln47r, as usual in the MS scheme. In eq. (|l^) the 
{A)Pij{z) denote the unpolarized (polarized) one-loop splitting functions for the 
transitions j -^ i |^. The functions {A)fij{z) represent the freedom in choosing 
a factorization prescription. In the MS scheme these functions vanish. Note that 
even in the polarized case only the unpolarized Hij contribute to the factorization 
of final-state singularities, as we do not consider the production of polarized 
hadrons. 

Before proceeding, we have to mention an important subtlety related to 
the use of the HVBM prescription for 75, which affects the polarized function 
AHgg. It is a well-known property of the HVBM-75 that it leads to helicity 
non-conservation at the qqg vertex in d dimensions, expressed by a non-vanishing 
difference of unpolarized and polarized d- dimensional LO quark-to-quark splitting 
functions, 

A disagreeable consequence of this is a non-zero first moment (x-integral) of 
the non-singlet NLO anomalous dimension for the evolution of polarized non- 
singlet quark densities, in obvious conflict with the conservation of the flavour 
non-singlet axial current [^, |19|, Q. At the same time, ( pO| ) is responsible for 
producing a result for the 0{as) correction to the Bj0rken sum rule [^ which 
disagrees with the one of |2^. These two effects turn out to be closely related, 
as they can be simultaneously removed by a factorization scheme transformation 



, 20|, generated by the term on the right-hand-side of (PQ). In other words. 



it is advisable, albeit not mandatory in a purely mathematical sense, to slightly 



deviate from the MS scheme in the polarized case by choosing (see also ||20| ) 

Afggiz) = -4Cp{l - z) (21) 
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in (0). The factorization scheme transformation defined by this equation has 
also been performed in the calculations of the spin-dependent NLO splitting func- 
tions [l^, |2^ and is thus respected by the available sets of spin-dependent NLO 
parton densities [1-4]. The "75-effect" described above has been known to occur 
in the HVBM scheme for quite some time [23, 2^, 0, 25 1 and is obviously a 



pure artefact of the regularization prescription chosen. Since furthermore physi- 
cal requirements such as the conservation of the non-singlet axial current serve to 
remove the effect in a straightforward and obvious way, results of NLO calcula- 



tions in "spin-physics" (like the ones of [|T9|, BO], or ours) are usually regarded as 



being "genuinely" in the conventional MS scheme only after this transformation 



has been carried out. The quantities A/^g, A/^^, and A/^^ in (|T^) will of course 
be set to zero, as in the usual MS scheme. Needless to say that in the unpolarized 
case (MS) one has fgg = fqg = fgq = fgg = 0. 

Another comment concerns the functions {A)Hgy needed for factorizing initial- 
state coUinear singularities from photon-splitting to a qq pair. As mentioned 
above, such singularities are absorbed into the "pointlike" part of the photon 
structure functions. Studies [|26|, ^ |2^ of the photon structure beyond LO have 
revealed that the MS-scheme photonic coefficient functions for the photon's DIS 
structure functions F^, gj exhibit a logarithmically singular behaviour at large 
X. Combining at NLO the "pointlike" parts of F^, gj with estimates for the 
"hadronic" component based on vector meson dominance (VMD) arguments, 
one encounters strongly negative results at large x, ruling out the use of intu- 
itive VMD ideas in the MS scheme. Instead, an appropriately adjusted ("fine 
tuned") non-VMD hadronic NLO input would be required in the MS scheme, 
substantially differing from the LO one, as the only means of avoiding unwanted 
and physically not acceptable perturbative instabilities for physical quantities 



like F2, gJ. In the unpolarized case the so-called DIS-y factorization scheme [26 
was introduced to avoid such "inconsistencies" . Here the idea was to absorb the 
photonic Wilson coefficient for F^ into the photon's quark densities by a factor- 
ization scheme transformation, hereby leaving the "hadronic" part untouched. In 
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[ p^ , this procedure was extended to the polarized case. It was found that after 
transforming to the DIS^ scheme, a pure VMD input can be successfully used for 
phenomenological analyses going beyond the LO. We will therefore specify the 
functions (A)/^^ to be used to transform to the DIS^ scheme. They read: 

1 — X 






Tr 
Tr 



[x' + (1 
{2x - 1) 



X] 



In 



X 



1 +6x(l 



X] 



In- 
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-1 +2(l-a;) 



(22) 



where Tr = 1/2. Of course, the choice of factorization scheme cannot affect the 
result for a physical quantity. In other words, in the unpolarized case, where all 
contributions can be consistently calculated to NLO, it does not matter eventually 
whether we use photonic parton densities defined in the DIS^ or the MS scheme, 
as long as we use NLO hard cross sections determined in the same scheme. In 
the polarized case however, we are not yet able to consistently include the NLO 
"resolved" contributions, as was pointed out several times before. Therefore, 
comparing the results for the direct part of the NLO cross section in the MS and 
the DIS-y schemes might indicate the uncertainty resulting from not performing 
a consistent NLO calculation. 

2.7 Final results 



For all processes the final partonic cross section can be cast into the form: 
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ln(l — V + vw) 

+ ci2^— - + cn 

I — w 



(23) 



Distributions in w like 5(1 — w), 1/(1 — w)+, etc. only occur for the subprocesses 
that are already present at the Born level. An expression similar to (|23| ) holds for 
the unpolarized case with, obviously, different coefficients Ci{v,w). We note that 
we have compared our unpolarized results to the ones presented in an analytical 
form in [§]. We found an almost complete overall agreement; however, there are a 
very small number of differences, some of which could be related to typographical 
mistakes. The only major discrepancies arise for the subprocesses 7g — > q{q'q') 
and 7g -^ q'iqq')- For the first, we believe that the result in [Q was accidentally 
presented in terms of the "crossed" process g7 -^ qiq'q')- For •yq — > q'{qq'), it 
seems that the result in [Q rather corresponds to 7g — > q'{qq')- Anyway, none 
of these small discrepancies turns out to have a significant numerical effect. The 
coefficients Cj(f , w) for the unpolarized and polarized cases are rather lengthy and 
will not be given here. They can be obtained in a Fortran code via electronic 
mail from Werner. Vogelsang@cern.ch. 

3 Numerical results 

Let us now present some first numerical results for the NLO corrections to po- 
larized single-inclusive photoproduction of charged hadrons. Rather than per- 
forming a detailed numerical study of the process, we will restrict ourselves to 
the most interesting questions. These concern the general size of the correc- 
tions ( "Ji'-f actor s" ) and the residual dependence of the NLO cross section on the 
unphysical scales present in the calculation. 

Before starting, we mention that whenever we will calculate the unpolarized 
NLO cross section, we will do so in a completely consistent way, i.e. by including 
both the direct and the resolved parts at NLO. Here we make use of our own 
results for the NLO corrections to the direct part of the cross section (see sec. 2), 
and of the ones in ||^ for the NLO resolved part. Furthermore, we will for 
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consistency use NLO parton densities for the incoming proton |2^ and the photon 
p^ , as well as NLO fragmentation functions. For the latter we will use the ones 
of [^ set up for the sum of charged pions and kaons. They will also be our choice 
when calculating the polarized cross section. 

In the polarized case at NLO, we will use spin-dependent parton distributions 
for the proton evolved at NLO and fitted to the available DIS data. Several sets 
for these are available |l|, ^ ^; for definiteness we will choose the ones of [|l| 
determined within the "radiative parton model". These also have the agreeable 
property of providing parametrizations at NLO and LO, the latter to be used 
for Born level predictions. In particular, we will choose the "valence" set of |[^, 
which corresponds to the best-fit result of that paper, along with one other set of 
[Q based on assuming Ag{x,fi'^) = g{x,^'^) at the low input scale /i of [Q, where 
gi^Xjfi"^) is the unpolarized GRV [^ input gluon distribution. This set will be 
referred to as "max. gluon" in what follows. Employing these two sets, which 
both provide a good fit to the available DIS data, but differ significantly in the 
polarized gluon density, we are able to see to which extent the relative size of the 
NLO corrections depends on the set of parton distributions used. 

We also note that whenever we calculate a cross section at LO (for instance, 
when determining the i^-factor K = a^^'^/o"^'^), we will for consistency use LO 
parton distributions and fragmentation functions. In this case we will also use 
the one-loop expression for the strong coupling, whereas at NLO we obviously 
employ its two-loop counterpart. The LO/NLO values for the QCD scale pa- 
rameter ^QCf) for Uf active flavours are taken from |^, 26, |l|. Heavy flavour 



(c, b) contributions to the cross sections are neglected for simplicity. Unless we 
exphcitly study the scale dependence of our results, we will choose the renormal- 
ization and factorization scales to be equal to the transverse momentum px of 
the produced hadron. 

We will provide numerical results for both the fixed target and the HERA 
collider kinematic domains. While the resolved component is expected to be 
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Figure 1: Ratio of LO direct and full (direct + resolved) polarized cross sections for 
HERA energies {Ef, = 27 GeV, Ep = 820 GeV). a. px dependence at 77 = —1, b. rj 
dependence at pr = 5 GeV. 

generally small at fixed target energies, it is known |^ to be dominant in certain 
regions of phase space at HERA also for the polarized case. Here the direct 
contribution will dominate only at fairly large pt, and/or at negative rapidities 
rj of the produced hadron in the HERA laboratory frame, where we have, as 
usual, counted positive rapidity in the proton forward direction. In order to 
demonstrate this, and to isolate for our further HERA studies the region where 
the direct contribution dominates, figure 1 shows the ratio of the direct part of 
the polarized cross section over the full (direct + resolved) one, calculated at LO 
and plotted vs. Pt {a^t rj = —1) and 77 (at Pt = 5 GeV). We have assumed E^ = 27 
GeV and Ep = 820 GeV; the cuts on the polarized Weizsacker- Williams spectrum 
were chosen as in |^. We have used the GRSV "max. gluon" set for the polarized 
proton. For the LO resolved part in the denominator we have to pick a suitable 
set of LO parton distributions for the polarized photon. Of course, nothing is 
known as yet experimentally about the latter, so we need to resort to models for 
them. Here we will follow ^ to use two very different scenarios, first considered 
in [RU|. They are based on assuming "maximal" (A/'^(x,/i^) = /'^(x, /i^)) or 
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"minimal" {Af^{x, ^'^) = 0) saturation of the fundamental positivity constraints 
|A/'^(x,/i^)| < p^x.jj?) at the input scale /x for the QCD evolution, where /x 
and the unpolarized photon structure functions f"^{x,fi'^) were adopted from the 
phenomenologically successful radiative parton model predictions in |^. These 
sets will be dubbed "max. 7" and "min. 7" sets, respectively, and figure 1 
shows the results obtained for both sets. As can be seen, in the region defined 
by r] < — 1, pt > 5 GeV the resolved component is expected to contribute about 
20% or less to the cross section (note that the direct and resolved parts of the 
cross section turn out to be of opposite sign). 

Having determined the region where the direct component dominates for 
HERA energies, we can now turn to NLO. Figure 2 shows the i^-factors for 
the direct part of the polarized cross section in the MS scheme, again vs. pt 
(at 7] = —1) and rj (at px = 5 GeV). The solid line corresponds to the "max. 
gluon" set for the polarized parton densities of the proton, whereas the dashed 
one displays the result obtained within the "valence" best-fit scenario of |l|]. As 
one can clearly see, the i^-factors are of very moderate size, K < 1 for almost 
all Pt and rj examined. Only at very large px, near the edge of phase space for 
the 7] = —1 considered, does the ii"-factor become much larger than unity within 
the "valence" scenario. This finding of generally small NLO corrections is very 
important and corroborates the LO predictions previously made in ||^. 

As frequently mentioned earlier, the NLO direct part on its own is factoriza- 
tion scheme dependent. For comparison we also plot in fig. 2 the /('-factor for 
the direct cross section obtained within the DIS^ scheme introduced in sec. 2.6. 
As can be seen, the corresponding change of the result is rather small. Finally, 
figure 2 also presents the if -factor for the full (direct -|- resolved) unpolarized 
cross section, which of course is scheme-independent. It turns out that it is very 
similar in size and shape to the K-factors we have obtained for the direct part 
of the polarized cross section. This, again, is a very satisfactory finding, as it 
suggests that our K-ia.ctoT for the direct part might not be too far off the result 
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Figure 2: ii'-factors for the direct part of the polarized cross section at HERA energies 
for different GRSV [^ parton distributions. In a. the K-factor for the total unpolarized 
cross section is also shown. 

for the one of the full polarized cross section, to be eventually determined when 
the NLO corrections to the resolved part of the polarized cross section will have 
been calculated. 

Another important issue when going beyond the LO is the expected reduction 
in the dependence of the results on the unphysical scales fi, M, Mp introduced 



previously. We now set /i^ = M^ 



Ml = ^pI and plot in fig. 3 the LO and NLO 



direct cross sections as functions of C, for fixed rj = —1, pt = 5 GeV. Even though 
we can only consider the direct part, the improvement in the scale dependence 
when going from LO to NLO becomes already clearly visible. 

We finally turn to the fixed target region, relevant for the HERMES and the 
future COMPASS experiments. It is again interesting to study the size of the K- 
factor for this situation, choosing a muon beam energy of 200 GeV. The results for 
our two sets of polarized parton densities of the proton are displayed as functions 
oi Pt in fig. 4, where we have fixed the centre-of-mass rapidity, rjcm = 0. We have 



19 



-10 - 



-11 




0.6 



1.2 



1.4 



1.6 



Figure 3: Scale dependence of the direct part of the polarized inclusive-hadron pho- 
toproduction cross section at p-p = 5 GeV and r^ = — 1 for HERA energies. All scales 
have been set equal to y/^pr, and the parton distributions used correspond to the 
"max. gluon" set of [p. The NLO cross section has been calculated in the MS 
scheme. 

again calculated the NLO cross section in the MS scheme. One can clearly see 
that again the if -factors are of very reasonable size, once Pt > 3 GeV, where one 
intuitively would start to trust perturbation theory. 

4 Summary and Conclusions 



We have presented for the first time the next-to-leading order QCD corrections 
to the spin-dependent cross section for single-inclusive charged-hadron photo- 
production. This process derives its importance from its sensitivity to the pro- 
ton's spin-dependent gluon distribution and, at high energies, to the so far com- 
pletely unknown parton content of circularly polarized quasi-real photons. It 
could be studied experimentally in future polarized fixed-target lepton-nucleon 
experiments, but also at the HERA ep collider after an upgrade to both beams 

20 



K 



1.5 



0.5 







llcm^O 


fixed target 


1 
1 








1 






GRSV max. gluon MS 


1 
1 
\ 







- - GRSV valence MS 


\ 

\ 
\ 




,---- 


- - ^ ^ 


\ 


^ -- 










^ 


^ ^^-<: 


r 

1 


^ 







PT[GeV] 

Figure 4: i^-factors for the direct part of the polarized single-inclusive charged-hadron 
cross section in a fixed target experiment with s = 400 GeV^ at r]cm = 0. 



being polarized. 

Our calculation is an important first step in trying to assess the perturbative 
stability of this process. First numerical results show generally moderate NLO 
corrections for the direct part of the cross section, the i^-factor being close to 
unity over a wide kinematical range at both HERA and fixed target energies. 
Also, the expected reduction in scale dependence of the cross section when going 
from LO to NLO is found. We finally emphasize, however, that in order to be able 
to use our results for obtaining truly physical predictions, the NLO corrections to 
the resolved part of the cross section will also have to be calculated in the future. 
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